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1.0 SUMMARY 


This doousant describes the Laval C requirements for the Shuttle Orbiter termi- 
nal area energy management (TAEM) guidance and flight control functions to be in- 
corporated into the Mission Control Center (MCC) entry profile planning 
processor. The logic presented is the required configuration for the first Shut- 
tle orbital flight (STS-1). TAEM requirements for flights subsequent to STS-1 
are presented in the optional TAEM targeting document (IN 80-FM-29). This 
processor will be used for preentry evaluation of the entry through landing 
maneuvers, and will lnolude a simplified 3-degree-of-freedom model of the body 
rotational dynamics that is necessary to aooount for the effects of attitude re- 
sponse on the trajeotory dynamics. This simulation terminates at TAEM-autoland 
interface. 

The TAEM guidance mode is initiat- based upon Em'** relative velocity. The 
TAEM guidance controls the Orbiter to altitude *-ic iyn*«ic pressure profiles 
that are Amotions of range -to-go. Angle of at. tact is r; 'ci \.u the alti- 

tude, and speedbrake deflection is used to control tiyr>e,v' preset tv.. Excess 
energy or insuffioient energy conditions result in u>,e a-.. turns cr o^layeu 
entry into the final approach phase, respectively. Roi- angle (a vaad for lat- 
eral trajectory control and for the S-turns, 

The body rotational dynamics simulation is a simp. 1 ' • cel t ,r.> assumes that 

the oommande.?" rotational rates oan be achieved within . - .nte^-ation interval. 
Thus, the rotational dynamios simulation is essentially a sl-'.uietion of the 
autopilot commanded rates and integration of these rates to determine the 
Orbiter attitude. 

It is assumed that certain data are available as input to this processor. It is 
assumed that the state vector is available in runway coordinates and the Orbiter 
heading is known with respect to the runway. Also, it is assumed that altitude 
of the vehicle center of mass with respeot to the runway, altitude rate, horizon- 
tal component of velocity, flightpath angle, and body roll and pitch attitudes 
are available as input data. Derivation of these data is presented in appendix 
A. 


2.0 INTRODUCTION 


This document defines the TAEM guidance and body rotational dynamics models 
required for the MCC simulation of the TAEM mission phase. This simulation 
begins at the end of the entry phase and terminates at TAEM-autoland interface. 

The TAEM guidance is simulated in detail. The rotational dynamics simulation is 
a simplified model that assumes that the commanded rotational rates can be 
achieved in the integration interval. Thus, the rotational dynamics simulation 
is essentially a simulation of the autopilot commanded rates and integration of 
these rates to determine Orbiter attitude. The rotational dynamics simulation 
also includes a simulation of the speedbrake deflection. The body flap and 
elevon deflections are computed in the Orbiter aerodynamic simulation. 
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3.0 OOOIPIHATB SYSTEMS 

The TAEM guidano* requires input data rafarenoad to two primary coordinate 
ayataaa. The first is tha topodatlo system whloh has its x-, y~, and z-axes 
oriantad in tha north, east, and down directions, respectively. Tha origin of 
this systaa is tha rehiole canter of aass, and the downward direction is noraal 
to the Barth ellipse id model. The seoond is the runway coordinate systaa. This 
systaa is obtained by constructing a topodatlo systaa with its origin at the oen- 
ter of the runway threshold and performing a rotation about its z-axls through 
the runway azimuth from true north. Both of these systens are Earth-fixed 
rotating systems, with the topodatlo being redefined at each time of interest. 

Pitch jnd roil attitudes are the angles of a yaw-pitoh-roll Euler rotation se- 
quence from topodatlo to body coordinates. 


4.0 TAEM GUIDANCE 

The TAEM fr.v.danoe will nominally guide the Orbiter to an altitude and dynamic 
pressure versus range profile. It oontrols energy by modulating drag through 
speedbrake deflection commands to null out the dynamic pressure errors and by 
nulling out altitude errors through the normal load factor oomraand. If excess 
energy exists, tha logic exeoutee an S-turn to dissipate additional energy. 

The previously available capability for selection of an alternate heading 
alinement cylinder nearer the runway (minimum entry point) haa been deleted and 
replaced by an option to force the selection of the left-hand heading alinement 
cylinder. This option is implemented by adding MED input capability for the 
existing minimum entry point flag (MEP) and by redefining its function as 
follows: 

MEP s 0, allow automatic selection by guidanoe of the nearest heading 

alinement cylinder (same side of runway oenterline as the Shut- 
tle position). 

MEP s 1 , force guidance to target for the left-hand heading alinement 
cylinder independent of Shuttle position. 

Note that this option does not provide selection of either heading alinement cyl- 
inder. For trajectories approaching from the left of runway centerline, both op- 
tions will result in a left-hand heading alinement turn. The left-hand HAC flag 
(MEP) will be initialized to zero in the mission constants table and changed by 
MED input if required. 

For all guidance phases (S-turn, acquisition, heading alinement, and prefinal), 
(IPHASE = 0, 1, 2, 3, respectively) the guidance normal load factor command is 
based on the altitude and altitude rate errors from the reference profile, and 
the speedbrake command is based chi the dynamic pressure error below >«'arh 0.9. 

All guidance phases are internally determined. 
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On fchs lateral axis, if th« guidanoe is in tha S~turn phase, s oenatant bank 
.angle is sent to the flight oontrol system. In the aoqulsltion phase, the guid- 
ance oo sound s a bank angle wiiioh is proportional to the Orbiter heading devia- 
tion from the tangent point on the nearest heading alinement oy Under. In the 
heading alineawnt phase, the guidance ooaoands a bank angle wlhioh assures that 
the Orbiter follows the heading allnement cylinder. During the prefinal phase, 
it ooomands a bank angle from a linear combinations of Orbiter lateral deviation 
and deviation rate from the runway centerline. 


4 . 1 DETAILED REQUIREMENTS 

The following subseotions 4.1.1 to 4.1,11 define the detailed software require- 
ments for the funotlons and subfunotlone that constitute the TAEM guidanoe pro- 
gram (ref, 1), These subseotions present the equations and the logio performed 
by function. The input and output variables are summarized in tables I and II. 
A summary of all constants is shown in table III. The values of the constants 
in this table are for a typieal OFT-1 trajectory, and are expected to change 
from mission to mission. 

The functions and sub functions of TAEM guidanoe are: 

TAEM guidance (TGEXEC) 

TAEM guidance initialization (TGINIT) 

XHAC function (TGXHAC) 

Groundtrack predictor (GTP) 

Resolve to 180 degrees routine (RES180) 

References and dynamic pressure function (TGCOMP) 

Phase transition and MEP function (TGTRAN) 

Normal acceleration command function (TGNZC) 

Speedbrake command function (TGSBC) 

Roll command function (TGPHIC) 

Functions and routines used in more than one other function or subfunction are: 
MID VALUE (MIDVAL) 


4.1.1 TAEM Guidance (TGEXEC) 

This function is the TAEM guidance executive routine. On the first pass, the 
initialization flag, IRESET, is set to 1, the initialization function (TGINIT) is 
executed first. Next the XHAC (TGXHAC) and the groundtrack predictor (GTP) func- 
tions are called. The TAEM reference and dynamic pressure function 'TGCOMP) is 
then executed. The TAEM phase transition and MEP function (TGTRANl is tn«n 
called. The guidance commands are then generated by calling the norcai accelera- 
tion command (TGNZC), the speedbrake command (TGSBC), and the roll command 
(TGPHIC) functions. TAEM guidance is then exited. 
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4 , 1 , 2 TABM Guidance Initialisation (TOIMITS 

The TOIMIT function initialises or computes uviml parameters used in TAEM guid- 
ance, Those operations are: 
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4.1.3 XHAC Function (TGXHAC) 

This function calculates the X runway coordinate of the heading alinement cylin- 
der center, XHAC, and the prefinal initiation range value (SHPLYK) as a function 
of surface wind conditions and weight. 

The operations performed are: 

XFTC = XA(IGI) + HFTCUGS)/TGGS(IGS) 

XALI = XA(IGI) + HALI ( I GS ) / TOGS ( I GS ) 

where the subscript IGI is selected based on a MED input GI_change and the sub- 
script IGS is automatically selected based on the input Orbiter mass. GI_change 
is a function of surface wind conditions. 

IGI s 1 if GI_change = 0 

IGI = 2 if GI change s 1 
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GX_ohang« Is sst to ztro if ths surface winds art loss than 50 percent of the de 
sign surfsos winds, end sst to 1 for winds greater than 50 psrosnt of ths design 
winds. 

Zf ths Orbiter aass is greater than WT_0S 1 , then I OS ■ 2 . Otherwise, IGS > 1. 
Ths value for WTJ3S1 is selected suoh that the glide slope ohange ooeurs when 
Orbitsr pay loadstars greater than 32 thousand pounds. 

Capability to seleot the minimum entry point (XALX) has been deleted, so XHAC is 
always set equal to XPTC. 


Ground track Computation (OTP? 

This section describes the equations used by the groundtraok predictor 
subprogram in mathematical symbols. The detailed formulation of the groundtraok 
predictor subprogram is given at the end of this section. 

An essential part of the TAEM guidance concept is the prediction of groundtraok 
range to runway thereshold. The method for computing groundtraok range is 
illustrated In figure 1. As shown, the Orbiter velooity vector is initially 
pointing away from tangency to the heading alinement cylinder. The acquisition 
turn is defined as the turn that will aline the Orbiter heading tangent to the 
heading alinement cylinder. The distance, d^Ci along the acquisition turn is 
computed from geometrical equations and an estimate of the acquisition turn ra- 
dius. The distance d^ is computed from geometrical equations relating the po- 
sition of the Orbiter after the acquisition turn to the tangent point on the cyl 
inder. The distance in the heading alinement turn is computed from the turn 
angle in the heading alinement phase and the alinement cylinder radius. The 
distance dj is a fixed distance from runway threshold of the alinement cylin- 
der . 

The groundtraok predictor estimates range-to-go by computing the segments 
^ACj d l > d HAC» 511(1 d T* These segments are added to give the tota) predicted 
range. The terms that no longer apply as the approach progresses are either set 
to zero or dropped from the summations. 

No additional range component is required for the S-turn phase because the 
predicted groundtraok during this phase is computed as if the S-turn were to 
cease immediately. 

During phase three (prefinal approach), the estimated range-to-go is based on a 
straight line to the runway threshold, thus 


RpRED = V* 2 7 y2 


The geometry for computing the segments d^c and di is given in figure 2. 
The angle is the difference between the Orbiter heading, \p , and the 

heading, \|>r, to the heading alinement circle tangency point. The vector from 
the Orbiter to the center of the heading alinement cylinder is computed as 
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*CIR ■ *HAC " * 

*CIR * YS0M R TURX * * 

The distanoe to ths csntsr of th# heading aline ment oirole la 
«CIR ■ ^CIR ♦ * Z CIR 

and the straight line distanoe to the tanganoy point is 

r tam « Yr 2 cir - R z turm 

which is then limited so that Rj^ > 0. 

The headingt 'i'C* t0 the oenter of the heading alineaent circle is given by 
4>C * tan” 1 (Tcir/Xcir) 

The heading to the tangency point is then 

v|>r = 4»C " YSGN tan’ 1 (RtURN/ r TAN) 
and the heading error is 




An acquisition turn radius is computed as a function of a predicted average bank 
angle during the turn. The turn radius, R^, and the error angle, A\b are 
used to compute the distance, d A C to be traveled during acquisition. The 
distance, d-j , to the tangency point from the end of acquisition is computed by 
using A^, Raci and R TA n* 

The radius of the acquisition turn, Rac» is approximated by 


RaC = V h V/(g tan <j> avg ) 


where = horizontal velocity. 
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♦xvo ■ predicted average bank angle. 

V * velooity Magnitude 
The acquisition turn aro is 
d JtC • *AC ! A *I 

Ths distance, dt is derived froa the geometry of figure 2. 

A » fyc <1. - cos tty)) 

® * &TAN “ SAC i aln (fy)l 
d! a VaTT^ 
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GROUHDTRACX PREDICTOR SUBPROGRAM 


If IPHASE « 3, 

RPRED « SORT (X 2 + Y 2 ) 

and then exit OTP. Otherwise, if IPHASE < 2, YSQM * SIGH (Y) (YSGH • + 1). The 
left hand HAC seleot flag is then tested, if MEP « 1, YSGH « -1. 

Now execute equation set 1: 

1.1 XCIR « XHAC -X 

1.2 YCIR * YSGN RTURN - Y 

1.3 RCIR * SORT (XCIR 2 + YCH 2 ) 

If IPHASE s 2, execute equation set 2: 

2.. PSHA = ARCTAN2 (XCIR, YSGH YCIR) RTD 

2.2 If PSHA < 0., PSHA s PSHA + 360 

2.3 RPRED s RCIR PSHA DTR - XHAC 

and then exit GTP. Otherwise (IPHASE = 0 or 1), oompute the straight line dis- 
tance to the HAC tangency point by 

IF: RCIR > RTURN 

THEN: RTAN = SORT (RCIR 2 - RTURN 2 ) 

ELSE: RTAN = 0 

The heading error (from a heading tangent to the HAC) is then computed by equa- 
tion set 3: 

3.1 PSC = ARCTAN2 (YCIR, XCIR) 

3.2 PST = (PSC - YSGN ARCTAN2( RTURN, RTAN)) RTD 
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3.3 PST • RES ISO (PST) 

3.4 DPSAC « RES ISO (PST- PSD 5 

The acquisition turn radius and aro length art than ooaputed in aquation sat 4: 

4.1 PHAYO > PHAVGC - PHAY03 MACH 

4.2 PHAVO « MIDYAL (PHAYO, PHAVOLL, PHAYOUL) 

4.3 RTAC ■ VH V/(0 TAM (PHAYO DTE)) 

4.4 ARCAC > RTAC ABS (DPSAC) DTR 

The range from the end of the aoquiaition turn to the HAC tangenoy point, RC, 
la next oomputed by equation set 5: 

5.1 A * RTAC (1 .-COS (DPSAC DTR)) 

5.2 B * RTAN - RTAC ABS (SIN (DPSAC DTR)) 

5.3 RC * SQRT (A 2 + B 2 ) 

The turn angle around the HAC is then defined by 
PSHA r ABS (PST) 

If (XCIR < 0 and ABS(Y) < 2. RTURN ) or 

(YSGN SIGN(Y) < 0.) , PSHA = 360 - PSHA 

The acquisition or S-turn predicted groundtrack is then computed by equation 
set 6: 


6.1 ARCHA * RTURN PSHA DTR 

6.2 RPRED s ARCAC + RC + ARCHA - XHAC 

GTP is then exited. 
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#,1,5 TASK Rsftrsnos and Dynamlo Pressure Function (TGCOMP) 

Th* TOCOMP function computes energy, altitude, fllghtpath angle, and dynamlo 
pressure references as a funotlon of the predicted range (APHID}. The energy 
reference (EM) consists of two linear segments which are functions of DtPRID. 
Energy slope and intercept for this linear funotlon are switched when DRPMSD Is 
less than the range switch point (BOW_SPT), The altitude reference (M ref) is 
a cubic curve which is tangent to the~eutoland steep glide slope at HP KID « 
»XALI , 

If HP HID is less than ~XALI, the altitude reference is defined by the autoland 
steep glide slope. If DHPHED la greater than PBRC, the altitude reference is a 
linear function of DHPHED, 

The reference fllghtpath angle is the slope of the altitude versus predicted 
range curve at the ourrent predicted range. The dynamic pressure reference is 
a two segment linear function of predicted range. 

This function also computes current energy, filtered dynamic pressure, filtered 
dynamlo pressure rate, dynamic pressure error, and the commanded equivalent 
airspeed. 


TGCOMP SUBPROGRAM 

Upon entering TGCOMP, the fm lowing computations are made: 

The current energy over weight is computed in equation set 1 : 

1.1 DRPRED « RPRED + XALI 

1.2 EOW * H + V V/(2G) 


The energy reference is computed in equations set 2. If DRPRED is less than EOW 
(IGS) , set I EL s 2. Otherwise, set I EL * 1, 

2.1 EN s EN_C1 (IGS, IEL ) + (DRPRED - RN1(IGS))EN_C2(IGS, IEL) 

If DPRED > PBRC(IGS), the altitude reference is computed with equation 4.1. 

4.1 HREF s PBHC(IGS) + PBGC (IGS) (DRPRED - PBRC (IGS)) 

If DRPRED < PBRC(IGS), the altitude reference is computed with equations 

4.2 and 4.3. 
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4.2 HREF » MALI (IQS) - TQGSUGS) DRPRED 

4.3 If { DEFIED >0.) HJtSF ■ HI IF * DEFIED DEFIED 

{OJB2C_C3ClM) + DEFIED 
CUSIC.C4(IGS)) 

Th® dynaalo praasur® raferenc# (QBREF) , altitud® ®rr*or, flightp^th angle 
rafaranoe, and ran®® amor from r®f#r®ne# altitude prof 11® are computed in 
aquation sat 5*. 

5.1 If (DRPRED < PBRCQ(IGS) ) , QBREF * HIDVAL (QBRUL(IGS) 

+ QBC2(IGS) DRPRED, 
QBRLL(IGS) , QBRUL(IGS) ) 

5.2 If (DRPRED > PBRCQ(IGS)), QBREF * MIDVAL (QBRLL(IGS) 

+ QBCKIGS) (DRPRED 
- PBRCQ(IGS)), 

QBRLL(IGS) , QBRML(IGS) ) 

5.3 HERROR = HREF - H 

5.4 If (DRPRED > PBRC(IGS)), DHDRRF = - PBGC(IGS) 

5.5 If (DRPRED < PBRC(IGS)), DHDRRF = - MIDVAL (-TGGS(IGS) 

+ DRPRED (2 CUBIC_ 

C3( IGS) + 3 CUBIC_ 

C4(IGS) DRPRED), 

PBGC (IGS), -TGGS(IGS) ) 

5.6 DELRNG * HERROR/DHDRRF 

Filtered dynamic pressure, filtered dynamic pressure rate, dynamic pressure 
error, and the commanded equivalent airspeed are computed in equation set 6: 
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6.1 QBARD > MIDVAL ( CQG* ( QBAR - QBARF) , - QBARDL, QBARDL) 

6.2 QBARF * QBARF + QBARD*DTG 

6 ’,3 OBD * CDBQD*QBD * CQDG*QBARD 

6.4 QBERR * QBREF - QBARF 

6.5 EAS CMD • 17. 1865*SQRT(QBREF) 


4.1.6 Phase Transition and MSP Funotlon (TOTRAH) 

This module determines all TAEM phase transition, and sets the flag TG_ END which 
will terminate the MCC simulation. Also, this funotlon ohecks for an S»turn 
(IPHASE a 0) situation. The Saturn initiation and termination is based on an 
energy error. 

The TAEM program will be entered in the acquisition phase (IPHASE * 1). In the 
normal situation, the Orb iter will remain in this phase until the heading 
alinement phase is initiated. However, if the energy state is too high, the 
S-turn phase will be initiated unless the predicted range is less than the 
minimum range allowed for an S-turn (RMINST) . The direction in which the 3-turn 
is made depends on which side of the runway centerline the Orblter is located 
and on the Orblter heading. The logic for selection of the turn direction is 
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S s ~1 left bank 
S = +1 right bank 


Termination of the S-tum occurs when the current energy falls below the refer- 
ence energy versus range line (EN). 

Transition to the heading alinement phase from the acquisition phase is tested 
whenever RCIR < P2TRNC1 RTURN where RCIR is the Orbiter's radial distance from 
the HAC enter and RTURN is the HAC radius. If this test is passed, transition 
to the heading alinement phase will occur if either of the following tests is 
passed: 
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a. IT < RTBIAS 

b. 1C IB < P2TRNC2 RTUIi 


where 


IT > (XDOT XCI1 * YDOT YCIJD/Y 

Transition to the prefinal phase CIPHASE * 3) oocurs whenever the predicted 
range is less than the prefinal Initiation range value (SHPLYK) . 

After phase 3 la initiated, the TAEM/ auto land transition ioglo is executed to 
determine when the transition to autoland will take place and the MCC simulation 
terminated . 

TGTRAH SUBPROGRAM 

If IPHASE = 3, a logical test is made for termination of TAEM guidance by the 
statement; 

If (IHERROR! < H_ERR0R and ( IY! < YJSRROR) and 

( ! GAMMA - GAMSGS(IGS) ! < GAMMAJSRROR) and 
( IQBERRI < QB_ERR0R 1 ) and 
(H > H_REF1 ) 

or 

( IHERROR I < (H DEL_H1 - DEL_H2 ) and 
( I Y l < (H Y_RANGE1 - Y_RANGE2) a..u 

(IGAMMA - GAMSGS ( IGS ) ! < (H GAMMA_C0EF1 - GAMMA_COEF2) and 
(IQBERRI < QB_J2RROR2) and 
(H_REF1 > H) 

or 


(H < H REF2) 
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then the TAEM guidanoe termination flag ia aet to one; i.e., TGJ5HD ■ 1, and 
then the TOTRAM function ia exited. " 


Otherwise , a teat for tranaition to the p re final phase ( I PHASE * 3) is made on 
the basis of whether the predioted range la leaa than the prefinal initiation 
range value (RPRED < SHPLYK). Zf thia teat la true, then equation aet 1 la 
exeouted: 


1.1 I PHASE * 3 

1.2 PMIO « PHIC 

1.3 PHILIM « PHILM3 
l.i) DNZUL « DMZUC2 
1.5 DNZLL * DNZLC2 


If the RPRED < SHPLYK teat ia falae, the appropriate logic baaed on the ourrent 
TAEM phaaa ia executed aa followa: 


For IPHASE * 0 (S-tum): 


If the current energy (EOW) ia greater than the reference energy (EH) exit 
TGTRAN . 

If the current energy ia leaa than the reference energy, then execute equation 
aet 2: 


2.1 

IPHASE s 1 

2.2 

PHILIM s PHILM1 

For IPHASE s 1 

(acquisition) 


Firat the energy over weight value to initiate an S-turn ia calculated aa shown 
in equation 3,1: 


3.1 ES = ESI(IGS) + ( DRPRED - RMINST(IGS)) EDRS(IGS) 


Next, the test for transition to the S-turn phase is made on the basis of pass 
ing the following criteria: 
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SOW > ES and DRPR2D > RMIM3T (IQS) 


If this tsst is passed, equation sst 3*2 is executed; 

3.2.1 I PHASE « 0 

3.2.2 PHILIM * PHILMO 

3.2.3 S « YSGH 

3.2.4 SPSI s S PSD 

3.2.5 If SPSI < 0 op SPSI > 135) then S » -S 

3.2.6 If MEP * 1, thsn S * 1. 

Next, ths energy over weight value used In the MEP teat is calculated as shown 
In equation 3.3: 

3.3 EMEP s EMEP_C1 (IGS, IEL) + (DHPRED - RH1(IGS))EHEP_C2(IGS, IEL) 

The test for transition to the heading alinement phase is made. If 

RCIR < P2TRNC1 RTURN , the parameter RT is calculated as shown in equation 3.4: 

3.4 RT s (XDOT XCIR + YDOT YCIR)/V 

and a test is made for passing either of the following conditions: 

3.4A RT < RTBIAS 

3.4B RCIR < P2TRHC2 RTURN 

If either 3.4A or 3.4B is satisfied, equation set 3.5 is executed. 

3.5.1 IPHASE = 2 

3.5.2 PHILIM = PHILM2 

For IPHASE = 2 (heading alinement): 

No calculations are presently required. 
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4,1.7 MopmI Aooslsratlon Co— and Punotlon (TOMZC) 

The function TOMZC oomputss the incremental nonut 1 acceleration oommand , NZC, to 
tha PCS for all TAEM phases . The oommand is generated from a computed altitude 
rate error term (HDERR), whioh is formulated from the altitude error input from 
TGCOMD?, a computed desired altitude rate term (HDREF), and the current, sensed 
altitude rate input from navigation. The upper and lower limits on the command 
are based on dynamio pressure and TAEM phase considerations. During an S-turn, 
a correction term, GCONT, to the basic turn coordination term oomputed in the 
PCS Is added to NZC. A filter limiting the rate of change of the unlimited com- 
mand (DMZC5 is Implemented to smooth the command for all TAEM phases except 
pre final. 

Upon entering TOMZC, the unlimited normal acceleration command (DNZC) and the 
upper and lower dynamic pressure limits QBMXNZ, QBMNNZ are computed using equa- 
tion set 1. 


1.1 GDH s MIDVAL (GDHC - GDHS H, GDHLL, GDHUL) 

1.2 HDREF s VH DHDRRF 

1.3 HDERR s HDREF - HDOT 

1.4 DNZC = DNZCG GDH (HDERR + HDREQC GDH HERROR) 

1.5 QBMNNZ = QBLL/AMAX1 (COSPHI, CPMIN) 

1.6 QBMXNZ * QBMX1 

1.7 IF MACH > QBM1 then 
QBMXNZ = MIDVAL ( QBMX2 

+ QBMXS(MACH - QBM2), QBMX2 , QBMX3 ) 


Next, the upper and lower limits based on the minimum and maximum dynamic pres- 
sure (QBNZUL, QBNZLL) are computed using equation set 2: 


2.1 QBNZUL a -(QBG1 (QBMNNZ - QBARF) - QBD) QBG2 

2.2 QBNZLL = - (QBG1 (QBMXNZ - QBARF) - QBD) QBG2 


If IPHASE a 3, then NZC = MIDVAL (DNZC, QBNZLL, QBNZUL) and then exit TGNZC. 
Otherwise, the filtered and limited normal acceleration command (NZC) is 
computed using equation set 3: 
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3.1 SKJUC » H + BMLMZ(XOS) MIDVAL ((DRPIBD 

- RH1 (IG3))/DELJUWAX<I0S), EDELC1, BDELC2) 

3.2 EMIN > EN - EDBLHZ(IGS) 

3.3 EQWfZUL > (GEUL GDHCEMAX . SOW) ♦ HDBIR)OEHDUL GDH 

3.4 EOWNZLL * (SELL GDH(EHIM - EON) * HDSRR)GEHDLL ODH 

3.5 DMZCD * (DNZC - DMZCF) COG 

3.6 DMZCF « MIDVAL (DMZCF ♦ DHZCD DTO, EOWNZLL, EOWMZUL) 

3.7 NZC « MIDVAL (DNZCP, QBNZLL, QBNZUL) 

* 

If I PHASE i 3, a test for the S-turn phase Is made (IPHASE * O); and if true, 
the turn coordination oorreotion ter® (GCONT) is oomputed and added to NZC using 
equation set 4: 

4.1 GCONT * (1. - QBARF/QBMXNZ) TAS/((TAS + VCO) SECTH COSPHI) 

4.2 NZC * NZC - QCONT 


Finally, for all phases, NZC is limited between the upper and lower phase depen- 
dent limits (DNZLL, DNZUL) by NZC * MIDVAL (NZC, DNZLL, DNZUL) . 


4.1.8 Speedbrake Command Function (TGSBC) 

This module performs the calculations required to generate the speedbrake com- 
mand to the FCS for all TAEM phases. The upper limit on the speedbrake command 
(DSBLIM) is Mach-dependent and is set to 65 degrees for supersonic flight and to 
98.6 degrees for subsonic flight. During the S-turn phase (IPHASE * 0), the 
speedbrake command is simply DSBLIM. For all other phases, the command is based 
on a nominal speedbrake command (DSBNOM) , plus proportional and integral terms 
based on the dynamic pressure error (QBERR) input from TGCOMP. The integral com- 
ponent of the command (DSBI) is held at last computed value whenever the previ- 
ous computed speedbrake command equals or exceeds either the upper or lower 
speedbrake command limit. 

If MACH is greater than DSBCM, the speedbrake command is set to a constant value 
of DSBSUP and TGSBC is exited. Otherwise, the TAEM phase is tested and if 
IPHASE * 0, equation set 1 is executed. 


1.1 DSBC s DSBLIM 
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Otherwise (IPHASE 4 0), aquation Mt 2 la executed. 


2.1 DSBB • QSBE QBE HR 

2.? If DSBC < DSBLIM and DSBC > 0.0, than 

D3BI ■ HIDVAL (DSBX * GSBI QBBNt DTQ.-DSBIL, DSBIL) 
2.3 DSBC > D3BNOM - DSBB - DSBX 


Finally, tha speedbrake oommand Issued to tha FCS, DSBC JIT, la oaloulatad using 
aquation aat 3; 


3.1 DSBC AT b MIDVAL (DSBC, SBHIM, DSBLIM) 


4,1.9 Boll Command Function (TGPHIC) 

Tha funotion TGPHIC oooputas the roll command, PHIC_AT, for the lateral axis com- 
asnd to the FCS for all TAEM phases. If the guidance is in the S-turn phase 
(IPHASE * 0), a oonstant roll command equal to S PHILIMIT is Issued to the FCS 
where S is the sign of the roll command calculated in TGTRAH , and PHILIMIT is 
the maximum roll command allowed. During the acquisition phase (IPHASE m 1), a 
roll command is given which is proportional to the Orbiter's heading deviation 
from tangsnoy to the HAC. In the heading alinement phase (IPHASE * 2), the roll 
oouffliand is generated to assure that tha Orbiter performs a turn which follows 
the heading alinement cylinder. In the prefinal phase (IPHASE * 3), the roll 
oowand is generated from a linear combination of the Orbiter's lateral devia- 
tion and deviation rate from the runway centerline. The roll command limit, 
PHILIMIT, is fixed for supersonic flight and is equal to the phase dependent 
roll command limit, PHILIM, for subsonic flight. 

Upon entering TGPHIC, a test on MACH is made for setting the roll command limit, 
PHILIMIT. For MACH > PHIM, PHILIMIT * PHILMSUP, and for MACH < PHIM, 

PHILIMIT s PHILIM. Next, the unlimited roll command, PHIC, is generated 
on the basis of the current TAEM phase as follows: 


1. IPHASE = 0 (S-turn): 

1.1 PHIC m S PHILIMIT 

2. IPHASE 2 1 (acquisition): 

2.1 PHIC s GPHI DPSAC 

3. IPHASE s 2 (heading alinement): 

3.1 RDQT g (XCIR XDOT + YCIR YD0T)/RCIR 
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3.2 RXRRC s MIDVAL (QKRCXR - RTURM) , -RRRRLM, RJRRLM) 

3.3 PHIC b ISOM (PKXP2C * RBRRC * ORDOT RDOT) 

4. I PHASE b 3 (pre final)* 

4.1 IERRC b MIDVAL (-01 I, -liRRLM, IIRflLM) 

4.2 PHIC b IE IRC - OX DO T YDOT 

4.3 Xf XSR ^ 0 than 

4.3.1 DPHI b (PHIC - PHI0)/I3R 

4.3.2 XSR b XSR - 1 

4.3.3 PHIC b PMIO ❖ DPMI 

4.3.4 PHIO b PHIC 


Finally, after the phase dependent unlimited roll command is oaloulated, the 
roll command to the PCS, PHIC_AT, is oaloulated uslns equation set 5: 

5.1 PHIC_AT b MIDVAL (PHIC, -PHILIMIT, PHILIMIT) 

4.1.10 Mid Value (MIDVAL) 

The MIDVAL funetion applies upper and lower limits to a variable. 


a. Detailed requirements.- The oall to function MIDVAL by the statement CALL 
MIDVAL (X, XLL , XUL) oommands the limiting operations: 


IF (X > XUL), THEN X * XUL 
IF (X < XLL), THEN X * XLL 

4.1.11 Resolve to 180 Degrees (RBS180) 

The function of RES180 is to resolve the input angle to lie within the range 
+ 180 degrees. 
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Calling RES 180 with tha statement 


AHQ • RES 180 UNO) 


oausas tha funotion RES 1 80 to axaouta tha loglo shown in flgura 3* 


4.2 SEQUENCING REQUIREMENTS 

TASK guidanoa ia to ba exeouted at a rata of 1.04 HZ (.96 aeo). Tha guidanoa 
da tall ad flow la shown in flgura 4, 


4.3 INTERFACE REQUIREMENTS 

TAEM guidanoa Input ano output parameters ara listed In table I and table II, re- 
spectively. 

It is assumed that service routines In the Mission Control Center will compute 
all of the required inputs to the TAEM guidanoa and flight control system 
simulations. Tha state vector in runway coordinates and the Orblter heading 
with respect to the runway oenterline are required input data. Also, required 
input data are the dynamlo pressure, true airspeed, altitude rate, geodetic alti- 
tude above the runway, horizontal component of ground relative velocity, body 
roll, and pitch angles. The dynamic pressure and true airspeed data require- 
ments will be satisfied by using the true airspeed and air density from the envi- 
ronment to compute dynamic pressure. Derivation of these data is presented in 
appendix A. 

It is also assumed that the output parameters consisting of predicted range, 

TAEM phase counter, equivalent airspeed command, altitude error from reference 
glide slope, filtered dynamic pressure, and the energy deficiency alert flag 
will be sent to the display processor. 

The output parameters consisting of the normal load factor command, roll angle 
command, and speedbrake angle command are sent to the TAEM digital autopilot. 


5.0 TAEM DIGITAL AUTOPILOT (TDAP) 


5.1 REQUIREMENTS OVERVIEW 

The TAEM digital autopilot (TDAP) is the first order simulation of the Shuttle 
Orbiter flight control system which performs the functions of controlling the 
body attitude and speedbrake responses during the TAEM guidance phase of the 
entry trajectory. Formulation requirements for this function are presented in 
figures 5 through 8, and are discussed in detail in the subsequent sections. 
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5.2 FUNCTION MODULES 

The module TDAP consists of tho following principal subfunction modules. 

a. Tho roll ohannol oaloulotas body roll rat# using tho roll attltudo and TAEM 
guidanoo roll anglo ooaaand. 

b. Tho pit oh ohannol oaloulataa body pltoh rata as a function of noraal load 
factor, attltudaa, truo alrapood, and guidanoo noraal load faotor ooaaand. 

o. Tho yaw ohannol oaloulatoa tho body yaw rots required to produoo ooordlnatod 
turns and null lateral load f so tors. 

d. Tho attitude integrator oonverts body rates to a rotation aatrlx and updates 
the body to Inertial attitude aatrlx. 

e. Tho speedbrake ohannol oaloulatoa speed broke actuator rates using the guid- 
ance deflection ooaaand and hinge moment constraints, and updates the 
speedbrake deflection. 

The following utility subroutines are used by 1'DAP: 


a. FILTER ...... generalized first order filter 

b. SMOOTH ...... guidanoo ooaaand smoother 

0. EIGEN ...... eigen vector rotation of an orthogonal matrix 

d. MIDVAL ...... applies upper and lower limits 

e. MAMP multiplies two matrices 

f. CROSS ...... vector oross product 

g. Unit ....... unit vector calculations 

h. AMIN1 select most negative argument 

1. AMAX1 . select most postive argument 

j. SIN ....... trigonometric sine 

k. COS trigonometric cosine 

l. MOD modular counter conversion 


Formulation is supplied in this document for the functions of TDAP and the first 
three utility subroutines in the above list. 
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5.3 INTERFACE REQUIREMENTS 

Inputs rsqulrsd by TDAP arm listed In tsbls IV. Outputs sre listed in tsbls V. 
Constants required are shown In tsbls VI. 


5. A SEQUENCING REQUIREMENTS 

TDAP is to bs executed st s rats of 2.06 HZ (.A8 sso) during the TAEM guidanos 
phaas aftsr guidanos sx scut ion. 


; v INITIAI I Z AT I ON REQUIREMENTS 

Ths internal variabls LOOP in TDAP oust ha vs the value zero on ths first pass ex- 
ecution . The spssdbrake deflection, body to inertial attitude mtrix, body roll 
rats, and yaw rats sust bs available at ths initial execution pass as ths last 
pass values from ths entry guidance phase so they nay be subsequently updated by 
TDAP. 


5.6 METHOD 

The MCC entry profile planning processor is a 3-degree-of-freedom simulation in 
the sense that it includes complete second order equations of motion for three 
dimensional translational components only. The rotational equations of motion 
are first order and assume that commanded body attitude rates are achieved 
instantaneously. This assumption eliminates the requirement to calculate 
rotational moments and accelerations for the purpose of integrating rates and at- 
titudes as would be done in a 6-degree-of-freedom simulation. The formulation 
for attitude control consists of calculating desired body rates using the guid- 
ance roll and load factor commands. This formulation is derived from the level 
C flight software requirements for the entry flight control system (reference) 
and provides virtually the same trajectory and maneuvering characteristics as a 
6-degree-of-freedom simulation. 

The first calculations upon entering TDAP are to obtain the trigonometric func- 
tions of pitch, roll, and angle-of-attack attitude angles which are used 
subsequently. 

The roll channel calculations are performed next by first calculating the 
scheduled roll rate gain and roll rate limit. The slow rate TAEM roll command 
is then smoothed by linear extrapolation to the faster TDAP rate. The roll atti- 
tude error is calculated and converted to a roll which is then limited. 

The pitch channel begins with calculation of the open loop coordinated turn 
pitch rate and smoothing it through a first order lag filter. The TAEM NZ com- 
mand is linearly smoothed, and the load factor bias for equilibrium turn compen- 
sation is calculated. The smoothed NZ command is converted to a C* command 
which will produce the pitch rate required to maintain a constant load factor at 
the desired level. The C* command and turn compensation terms are summed to pro- 
duce the total desired load factor, and load factor error is calculated. 
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The load factor error is oonvtrted to a negative plfcoh rate command, smoothed 
through a first order lag filter, sunned with the coordinated turn pltoh rate, 
and filtered again to obtain the unlimited negative pltoh rate. Angle-of-attaok 
Units are obtalnsd from their scheduled profiles, and the maximum allowable 
pltoh rates are determined by the proximity from the angle~of*»attaok limits and 
load faotor limits. The aoat restrictive of these Halts are then applied to ob- 
tain the final pltoh rate. 

The yaw channel begins with calculation of the lateral load faotor yaw rate 
gain. The yaw rate required to maintain a ooordlnated turn maneuver about the 
X- stability ax' . *s determined from body roll rate and angle of attaok. The 
sensed body lateral load faotor Is smoothed through a first order filter, 
converted to a yaw rate, and summed with the coordination term to obtain the 
final yaw rate. 

To integrate the vehiole attitude, body rates are assumed constant over the com- 
putation interval and, therefore, form the eigen veotor in body coordinates 
about whloh the body to inertial attitude matrix is to be rotated. This matrix 
rotation is performed in the subfunotlons EIGEN by constructing the transforma- 
tion matrix whioh, when premultiplied by the present attitude matrix, will pro- 
duce the direction cosine matrix that represents the attitude resulting from an 
eigen rotation througi the total angle determined by the desired rate times the 
time interval. The columns of this matrix, which represent the updated body 
axes in the inertial frame, are then normalised to preserve orthogonality. 

The final function of TDAP is the first order simulation of speedbrake actuator 
response to the guidance commanded deflection. The principal considerations in 
determining the actuator rate include the flight software priority rate limiting 
and actuator ability to overcome the aerodynamic hinge moment. Since the 
Orbiter speedbrake actuator system gain is sufficient to produce time constants 
less than the TDAP cycle interval, a gain of 1/DT2 is used to prevent an 
oscillatory characteristic. The maximum opening rate is determined as a square 
root function of hinge moment. Maximum closing rate is constant at a value 
seleoted by the defleotion in relation to the software soft stop position. The 
deflection error is converted to a rate and limited. This rate is integrated 
rectangularly to obtain the speedbrake deflection which is then limited to the 
position constraints. 

The first order filter subfuncti^i utilizes the nodal implementation of the 
generalized S-transform (C-jS + CgJ/CS + a). Conversion to the equivalent Z~ 
transform is performed on the initial pass to obtain the required coefficients 
for the recursive difference equation. 

Detailed formulation for the function modules TDAP, FILTER, SMOOTH, and EIGEN is 
shown in figures 5 through 8. 
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6.0 R8FSI1WC1 

1. Spaoa Shuttla Orbital Flight Taat Laval C Functional Sub ays baa Softwara 
laqulrananta , 

Part A - Qui Sanaa, SO 76-3H-OO01C, Dao. 15, 1978. 

Part C ■ Flight Control - Entry, SO 76«SH«0007, Mov. 26, 1978. 
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TABLE I.- TAEM GUIDANCE INPUTS 


Description 

Type 

Souroe 

Unit 

Internal 

name 

Geodetic altitude of vehlole 
center of mass above RW 

F 

GCOMP* 

ft 

H 

Negative z component of velocity 
In topodetlo coordinates 

F 

GCOMP 

fps 

HDOT 

x~oomponant of position in runway 
coordinates 

F 

GCOMP 

ft 

X 

y»co«ponent of position in runway 
coordinates 

F 

GCOMP 

ft 

Y 

Magnitude of Earth relative 
velocity vector 

F 

GCOMP 

fps 

V 

Topodstio horizontal component 
of Earth relative velocity 

F 

GCOMP 

fps 

VH 

x-component of velocity in runway 
coordinates 

F 

GCOMP 

fps 

XDC T 

Course WRT RW centerline 

F 

GCOMP 

deg 

PSD 

MACH number 

F 

GCOMP 

— 

MACH 

Dynamic pressure 

F 

GCOMP 

psf 

QBAR 

Cosine of body roll Euler angle 

F 

GCOMP 

— 

COSPHI 

Secant of body pitch Euler angle 

F 

GCOMP 

— 

SECTH 

Orblter mass 

F 

AERO 

ibm 

WEIGHT 

True airspeed 

F 

GCOMP 

fps 

TAS 

y-component of velocity in 
runway coordinates 

F 

GCOMP 

fps 

YDOT 

a Guidance and flight control input parameter 

calculations 

(appendix) 
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TABLE I.- Concluded 


Description 

Type 

Souroe 

Unit 

Internal 

name 

Topodetio Earth relative flight* 
path angle 

F 

GCOMP 

deg 

GAMMA 

Flag indloating glide slope 
desired 

I 

HED b 

SMI 

GI_CHANGE 

Left-hand MAC select flag 

I 

MED 0 

* 

HEP 

b MED input based on surfaoe wind 

conditions 

(initialized 

in mission 

constants 


table to zero). 

°MED input equal to 1 to force selection of the left-hand heading alinement 
cylinder (initialized in mission constants table to zero). 
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TABLE II,- TAEM ©UIDAMCI OUTPUTS 


Description 

Type 

Unit 

Internal 

name 

Destination 

Commanded body normal load factor 
Increment from equilibrium 

P 

g's 

NZC 

TDAP 1 

Roll angle ooaaand 

F 

deg 

PHIGJIT 

TDAP 

Range error from referenoe 
altitude profile 

F 

ft 

DBLRHO 

DSP b 

Heading error to HAC tangenoy 

F 

deg 

DPSAC 

DSP 

Speedbrake angle command 
(hinge line) 

F 

deg 

DSBC_AT 

TDAP 

E/W at which the MEP is selected 

F 

ft 

EMEP 

DSP 

Energy per unit weight 

F 

ft 

EOW 

DSP 

E/W at which an S-turn is 
initiated 

F 

ft 

ES 

DSP 

Predicted range to runway 
threshold 

F 

ft 

RPRED 

DSP 

TAEM phase counter 

I 

— 

IP EASE 

DSP 

TAEM guidance termination flag 

I 

— 

TG_END 

DSP C 

Equivalent airspeed command 

F 

knots 

EAS_CMD 

DSP 

Altitude error from reference 
profile 

F 

ft 

HERROR 

DSP 

Filtered dynamic pressure 

F 

psf 

QBARF 

DSP 


a TAEM digital autopilot, 

^Displays. 

c This flag will terminate the MCC simulation. 
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TABLS III.- TAEH QUID AH CE CONST AMTS 



Description 

Value 

Unit 

CDSQD 

Gain in calculation of QBD 

.61113130 


CFMIH 

Minimum value CosFHI 

.TOT 


CQDG 

Gain used in calculation of 
QBD 

« 31 8S6S6Q 


CQG 

Gain used in oaloulstion of 
DMZCD end QBARD 

.5583958 


CUBIC — C3(l )* 

Quadratic tern eoaffletent 
in HREF CURVE FIT 

-.3641168E-6 

ft/ ft 2 

CUBIC_C3£2) 

Quadratic term coefficient 
in HREF CURVE FIT 

-.36411681-6 

ft/ ft 2 

CUBICj:4(1) 

Cubic term coefficient in 
HREF CURVE FIT 

-.9481026E-13 

ft/ft3 

CUBIC_C4(2) 

Cubic term coefficient in 
HREF CURVE FIT 

-.9481026E-13 

ft/ft3 

DELPHI 

Altitude error coefficient 

0.19 

— 

DELJJ2 

Altitude error coefficient 

900. 

ft 

DEL_R_EMAX( 1 ) 

Constant used in computing 
E max 

54000. 

ft 

DEL_RJMAX(2) 

Constant used in computing 
E max 

54000. 

ft 

DUZCG 

Gain used to compute DNZC 

0.01 

g/s 

DNZLC1 

Phases 0, 1, and 2 NZC 
lower limit 

-0.5 

g 

DNZLC2 

Phase 3, NZC lower limit 

-0.T5 

g 

DNZUC1 

Phase 0, 1, and 2 NZC 
upper limit 

0.5 

g 


a Array dimension 2. 
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TABLE ZIl.o Continued 


Symbol 

Description 

Value 

Unit 

DNZUC2 

Phase 3 NZC upper limit 

1.5 

8 

DSBCM 

Mach at which speedbrake aodulation 

begins 

0.9 


DSBLIM 

Maximum value of DSBC 

98.6 

deg 

DSBSUP 

Supersonic fixed speedbrake 
deflection 

65. 

deg 

DSBIL 

Limit on DSBI 

20. 

deg 

DSBNOM 

Nominal speedbrake command 

65. 

deg 

DSHPLY 

Delta range value in SHPLYK 

4000 . 

ft 

DTG 

TAEM guidance cycle interval 

0.96 

sec 

DTE 

Degrees -to- radian conversion factor 

0.0174533 

rad/deg 

EDELC1 

Constant used in determination of 
EMAX 

1 . 

HD 

EDELC2 

Constant used in determination of 
EMAX 

1 . 

HD 

EDELNZ(I) 

Energy delta from the nominal energy 
line slope for the S-turn 

4000. 

ft 

EDELNZ(2) 

Energy delta from the nominal energy 
line slope for the S-turn 

4000. 

ft 

EDRS(1) 

Slope for S-turn energy line 

.6089492 

ft 2 /ft 

EDRSC2) 

Slope for S-turn energy line 

.6089492 

ft 2 / ft 

EMEP_C1 ( 1 , 1 ) b 

Y-intercept of MEP energy line 

-271.3797 

1 „ 

EMEP_C1 (2 , 1 ) 

Y-intercept of MEP energy line 

-271.3797 

ft 

EMEP_C1 (1,2) 

Y-lntercept of MEP energy line 

1^604.87 

ft 


^2x2 array given column wise. 
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TABLE 111,- Continued 


Hi 


Value 

Unit 

EMEP_C1(2,2) 

Y-int«ro«pt of, MSP mwgy 11a* 

14604,87 

ft 

MPJC2(1,1) 

Slop* of MSP energy 11a* 

.4161274 

ft 2 / ft 

EHEP_ C2(2,1) 

Slop# of MEP energy line 

.4168274 

ft 2 / ft 

BCP_C2(1 t 2) 

Slop* of MEP energy lln* 

.2821187 

ft 2 / ft 

EMEP_C2(2,2) 

Y-intereept of nominal energy lln* 

.2821187 

ft 2 /ft 

EM_C1 (1,1) 

Y~interoept of nominal energy lln* 

-254.2395 

ft 2 /ft 

EN_C1(2,1) 

Y-lnt«ro*pt of nominal energy lln* 

-254.2395 

ft 2 /ft 

EN_C1 (1,2) 

Y-lntero*pt of nominal energy lln# 

17645.34 

ft 2 / ft 

EN _C1(2,2) 

Y-infc«rc*pt of nominal energy lln* 

17645.34 

ft 2 /ft 

BN_C2 (1,1) 

Slop* of nominal energy lln* 

.5500275 

ft 2 / ft 

EH_C2(2, 1 ) 

Slope of nominal energy line 

.5500275 

ft 2 /ft 

EH_C2(1,2) 

Slope of nominal energy line 

.3890240 

ft 2 /ft 

EN_C2(2,2) 

Slope of nominal energy line 

.3890240 

ft 2 /ft 

ES 1(1) 

Y-intercept of S-tum energy line 

90000. 

ft 

ES1(2) 

Y-intercept of S-turn energy line 

90000. 

ft 

E0W_SPT ( 1 ) 

Range at which to change slope and 
intercept on the MEP and nominal 
energy line 

117718. 

ft 

EOW_SPT ( 2 ) 

Range at which to change slope and 
intercept on the MEP and nominal 
energy line 

117718. 

ft 

G 

Earth gravitational constant 

32.174 

ft/s 2 

GAMMA_COEF 1 

Fllghtpath error coefficient 

0.0007 

deg/ft 

GAMMA_COEF2 

Flight path error coefficient 

3.0 

deg 
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TABLE III.* Continued 


Symbol 

Description 

Value 

Unit 

gmwajeiior 

Flightpath error band 

4.0 

deg 

OAMSGS(I) 

A/L eteep gli deslope angle 

-20 « 

deg 

OAHSOSC2) 

A/L eteep glidaslope angle 

-20, 

deg 

QDKC 

Constant for computing GDH 

2.0 

— 

GDHLL 

GDH lower Unit 

0.3 


GDHS 

Slope for computing GDH 

7.01-5 

ft" 1 

GDHUL 

GDH upper limit 

1.0 

— 

GEHDLL 

Gain used in ooaputing EOWNZLL 

0.01 

g/fps 

GEHDUL 

Gain used in computing EOWNZUL 

0.01 

g/fps 

CELL 

Gain used in ooaputing EOWNZLL 

0,1 

sec -1 

GEUL 

Gain used in ooaputing EOWNZUL 

0.1 

sec" 1 

GPHI 

Heading error gain for computing 
PHIC 

2.5 

— 

GR 

Gain on RCIR In computing HA roll 
angle command 

0.02 

deg/ ft 

GRPQT 

Gain on dRCIR/dt in computing HA 
roll angle command 

0.2 

deg/fps 

GSBE 

Speedbrak® proportional gain on 
QBERR 

1.5 

deg/psf 

GSBI 

Gain on QBERR integral in computing 
epeedbrake command 

0.1 

deg/psf-s 

GY 

Gain on Y in computing PPL roll 
angle command 

0.05 

deg/ ft 

GYDOT 

Gain on YDOT in computing PPL roll 
angle command 

0.6 

deg/fps 

H_ERR0R 

Altitude error bound 

1000. 

ft 
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Spfeol 

Description 

Value 

Unit 

xjnri 

Altitude reference 

2000 , 

ft 

H_§EP2 

Altitude reference 

2000. 

ft 

mid) 

Altitude et A/L eteep glidselops et 
MSP 

1001 a. 

ft 

HALI(2) 

Altitude et A/L steep glideslope et 
HEP 

10018. 

ft 

HDREQG 

Sain used to compute HDREQ 

0.1 

**>«-* 

HFTC(I) 

Altitude of A/L steep glideslope et 
nominal entry point 

12018. 

ft 

HFTCC2) 

Altitude of A/L steep glideslope et 
nominal entry point 

12018. 

ft 

MACHAD a 

Maoh no. to use elr date 

2.5 


MXQBWT 

Max L/D dynamic pressure for nominal 
weight 

0.73684218-3 

psf/lb 

PBGC CD 

Linear coefficient of range for HREF 
and lower limit for DHDRRF 

.1125953 

— 

PBGC (2) 

Linear coefficient of range for HREF 
and lower limit for DHDRRF 

.1125953 

— 

PBHCC1) 

Altitude reference for DRPRED * PBRC 

84821.29 

ft 

PBHCC2) 

Altitude reference for DRPRED * PBRC 

84821.29 

ft 

PBRC(1) 

Maximum range for cubic altitude 
reference 

308109.5 

ft 

PBRCC2) 

Maxim.. . range for cubic altitude 
reference 

308109.5 

ft 

PBRCQ(I) 

Range breakpoint for QBREF 

122000. 

ft 

PBRCQ(2) 

Range breakpoint for QBREF 

122000. 

ft 


a Usad to GCOMP calculations (appendix). 
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TABLE III.- Continued 


Symbol 

Description 

Value 

Unit 

P HAVOC 

Constant for computing PHAVO 

63.33 

deg 

PMAVOLL 

Lower limit for PHAVO 

30. 

deg 

PHAVOS 

Slops for oompucing PHAVO 

13.33 

deg 

PHAVOUL 

Uppsr limit for PHAVO 

50. 

deg 

PHILMO 

Roll oommand limit Saturn 

50. 

deg 

PHXLM1 

Roll command limit acquisition 

phass 

50. 

deg 

PHILM2 

Roll command limit heading 
alinement phase 

60. 

deg 

PHILM3 

Roll command limit pre final phase 

30. 

deg 

PHILMSUP 

Supersonic roll angle command limit 

30. 

deg 

PHIM 

Mach at which supersonic roll 
command limit is removed 

1 . 

— 

PHIP2C 

Constant for computing phase 2 roll 
command 

30. 

deg 

P2TRNC1 

Constant used in Phase 2 transition 
logic 

1.1 

— 

P2TRNC2 

Constant used in Phase 2 transition 
logic 

1.01 

— 

QB_ERR0R1 

Dynamic pressure error bound 

- 1 . 

psf 

QB_ERROR2 

Dynamic pressure error bound 

24. 

psf 

QBARDL 

Limit on QBARD 

20. 

psf/s 

QBCK1) 

Slope of QBREF for DRPRED > PBRCQ 

.5869565E-3 

psf/ ft 

QBC1(2) 

Slope of QBREF for DRPRED > PBRCQ 

. 5869565E-3 

psf/ ft 

QBC2C1) 

Slope of QBREF for DRPRED < PBRCQ 

~. 1031695E-2 

psf/ ft 
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TABLS III.- Continued 


Symbol 

Description 

Value 

Unit 

QBC2(2) 

Slop* of QBREF for DRPRED < PBRCQ 

-.1031695E-2 

pef/ft 

QB01 

Gain used to ooapute QBNZUL and 
QBNZLL 

0.1 

s" 1 

QB02 

Gain us ad to ooapute QBNZUL and 
QBNZLL 

0.125 

s-g/psJ 

QBHX3 

Slope of QBMXNZ with MACH > QBM2 

20. 

psf 

QBMX1 

Constant for ooaputing QBMXHZ 

340. 

paf 

QBMX2 

Constant for ooaputing QBMXNZ 

220. 

psf 

QBMX3 

Constant for ooaputing QBMXNZ 

250. 

psf 

QBM1 

Mach breakpoint for ooaputing 
QBMXNZ 

1.0 

e»<-» 

QBM2 

Mach breakpoint for computing 
QBMXNZ 

1 .0 

*a»*+ 

QBRLL(I) 

QBREF lower limit 

153. 

psf 

QBRLL(2) 

QBREF lower limit 

153. 

psf 

QBRML(I) 

QBREF middle limit 

180. 

r>sf 

QBRHLC2) 

QBREF middle limit 

180. 

psf 

QBRUL(l) 

QBREF upper limit 

265.43 

psf 

QBRULC2) 

QBREF upper limit 

265.43 

psf 

RERRLM 

Limit of RERRC 

50. 

deg 

RFTC 

Roll fader time constant 

5. 

S 

RMXNST(I) 

Minimum range allowed to initiate 
S-turn phase 

122204.6 

ft 

RHINST(2) 

Minimum range allowed to initiate 
S-turn phase 

122204.6 

ft 
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TA8LB III,- Concluded 


Symbol 

Description 

Value 

Unit 

RMU15 

itange at whioh the fain on XD3ELMZ 
sou to zsno 

65*2.886 

ft 

RHK2) 

Rang* at which tha sain on BDELWZ 
goes to ZIBQ 

65*2.885 

ft 

RTBIAS 

Maximal value for RT for HA 
initiation 

3000 . 

ft 

RTD 

Radians-to-degrees oonvaraion faotor 

57 .2957$ 

dag/rad 

RTU1M 

HAG radius 

20000. 

ft 

SBMIH 

Minimum speedbrake command 

5. 

deg 

TOOS(I) 

Tan of stesp gilds slops for autoland 

-.363978 

— 

TOGS (2) 

Ten or* stssp gilds slops for autoland 

-.363978 


VCO 

Constant used in turn compensation 

1.0120 

fps 

WTJ3S1 

Maximum Shuttle Orbitsr weight 

250000. 

lb 

XA( 1) 

Stssp glide si ops lntsrospt 

-6500, 

ft 

XA(2) 

Stssp glidsslops lntsrospt 

—5500 • 

ft 

YERRLM 

Limit on YERRC 

120. 

deg 

Y_ERR0R 

Crossrange error bound 

1000. 

ft 

YJtAHGEI 

Crossrangs oosfflolent 

0.18 

— 

Y RMGE2 

Crossrange ooeffloisnt 

800. 

ft 
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TABLE XV.- TDAP INPUTS 


Him 

Description 

Source 

Unit 

THETA 

Topodetio to body pitch angle 

OCOMpk 

deg 

PHI 

Topodetio to body roll angle 

OCOKP 

deg 

NZC 

Commanded body normal load faotor 
increment from equilibrium 

TAEM 

8 

TAS 

True airspeed 

GCOMP 

fps 

HZ 

Body normal load faotor 

GCOMP 

8 

MACH 

•iach number 

GCOMP 


ALPHA 

Angle of attack 

GCOMP 

deg 

NY 

Body lateral load faotor 

GCOMP 

8 

PHIC 

Roll angle command 

TAEM 

deg 

DSBC 

Speedbrake command (hinge line) 

TAEM 

deg 

DELB a 

Speedbrake deflection 

AERO 

deg 

HMSB 

Speedbrake hinge moment 

AERO 

in-lb 

CBI a 

Body to inertial attitude matrix 

EDAP C 

— 

P a 

Body roll rate 

EDAP 

deg/sec 

R a 

Body yaw rate 

EDAP 

deg/sec 


a Required only for initialization. 

^Guidance and flight control input parameter calculations. 
c Entry guidance phase autopilot. 
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TABUS Y.« TDAP OUTPUTS 


Mum 

Description 

Destination 

Unit 

P 

Body roll rite 

MOM* 

deg/seo 

Q 

Body pitoh rite 

MOM 

deg/seo 

R 

Body yaw rite 

MOM 

deg/seo 

CBI 

Body to inertial attitude matrix 

GCOMP 


DELB 

Speedbrake deflection (hinge line) 

AERO 

deg 

DDOTRC 

Speedbrake rate 

MOM 

deg/sec 


•■Not presently used except for possible oheokout monitoring. 
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TABLE VI.- TDAP CONSTANTS 



Description 

Value 

Unit 

AMNI 

T-interoept of minimum alpha profile 

**16, 

deg 

AMMLL 

Minimum alpha lower limit 

-100. 

deg 

hMS 

Slope of minimum alpha profile 

20. 

deg.' 

AHNUL 

Minimum alpha upper limit 

-100. 

dfc t 

AHXLL 

Maximum alpha lower limit 

100. 

deg 

AMXtft, 

Maximum alpha middle limit 

100. 

deg 

AMXM1 

Maximum alpha profile first 
breakpoint 

0.6 

mach 

AMXM2 

Maximum alpha profile second 
breakpoint 

0 . 

mach 

AMXS1 

Lower slope of maximum alpha profile 

-25. 

deg/mach 

AMXS2 

Upper slope of maximum alpha profile 

10. 

deg/ mach 

AMXUL 

Maximum alpha upper limit 

100. 

deg 

CPMIN 

Lower limit on cosine PHI 

0.5 

— 

DSBLIM 

Maximum speedbrake deflection 

98.6 

deg 

DT2 

TDAP time interval 

0.48 

sec 

DTG 

TAEM guidance interval 

0.96 

sec 

DTR 

Degrees to radians conversion 

0.0174533 

rad/deg 

GBDOT 

Gain to scale compensated yaw rate 

1 . 

— 

GF1 1 

Filter 1 first coefficient 

0 . 

— 

GF12 

Filter 1 second coefficient 

1 . 

— 

GF13 

Filter 1 third coefficient 

1 . 

— 

GF21 

Filter 2 first coefficient 

0.2 
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TABLE VI.- Continued 


Hants 

Description 

Value 

Unit 

GF22 

Filter 2 second coefficient 

2. 


GF23 

Filter 2 third ooeffloient 

2. 

~ 

GF31 

Filter 3 first coefficient 

2. 


GF32 

Filter 3 second coefficient 

2. 

— 

GF33 

Filter 3 third coefficient 

2. 

MM 

GF41 

Filter 4 first coefficient 

0. 


GF42 

Filter 4 second coefficient 

5. 


GF43 

Filter 4 third coefficient 

5. 

MM 

GPI 

Y-intercept of GPBANK profile 

4.4 

(deg/s)/deg 

GPLL 

GPBANK lower limit 

0.5 

(deg/s)/deg 

GPS 

Slope of GPBANK profile 

-3.25 

( deg/ s ) / ( deg-mach ) 

GPUL 

GPBANK upper limit 

1.8 

(deg/s)/deg 

GQA 

Alpha limit pitch rate gain 

1 . 

(deg/s)/deg 

GQN 

Pitch rate gain 

3-36 

(deg/s)/g 

GQNL 

Load factor limit pitch gain 

4. 

(deg/s)/g 

GRPHI 

Cross-over velocity gain 

-1845.06 

(deg/ft)/g-sec^ 

GY I 

Y-intercept of GRAY profile 

23.0 

(deg/s)/g 

GILL 

GRAY lower limit 

3.2 

(deg/s)/g 

GYS 

Slope of GRAY profile 

- 6.8 

( deg/s) /(g-mach) 

GYUL 

GRAY upper limit 

10. 

( deg/s )/g 

HMS 

Speedbrake stall hinge moment 

2.40E6 

in-lb 

NZMAX 

Maximum normal load factor 

100. 

g 
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TABLE VI t m Concluded 



Description 

Value 

Unit 

NZMIN 

Minimum normal load factor 

-100. 

S 

PCI 

Y«lnteroept of PCLIM profila 

30. 

deg/seo 

POLL 

PCLIM lower limit 

5. 

deg/sec 

PCS 

Slops of PCLIM profila 

-16.667 

deg/(sec-maon) 

PCUL 

PCLIM upper Halt 

20, 

deg/soc 

SBMIN 

Minimum epeedbrake deflection 

5. 

dag 

SBRLC 

Speedbrake maximum closing rate 

10.86 

deg/sec 

SBRLO 

Speedbrake maximum opening rate 

6.1 

dag/sec 

SBRS 

Speedbrake maximum closing rate 
below soft stop 

1 . 

deg/sec 

SBSOFT 

Speedbrake soft stop 

12. 

deg 

TPLIM 

Upper limit on tangent PHI 

1 . 

— 

VCO 

Cross-over velocity 

549.125 

fps 
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Runway 

threshold 



(Computed groundtrack range) 

Figure 1.- TAEM guidance groundtrack predictor geometry. 
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ENTER 

RES 180 

! 

i 

i 




\ 

i 

DO 

WHILE 

\ 

i ANG = ANG - 360. 

AND 

> 180. 

r 

1 

1 


"T 

i 

i 

t 

i 

i 

i 

i 

/ 




\ 

i ; 

DO 

WHILE 

\ 

I ANG = ANG + 360. \ 

AND 

< - 180 

/ 

/ 

I . 

1 J 


EyiT 


Figure 3.- Resolve to 180 degrees (RES 180). 
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TGEXEC 


8CTW25 


* 

t 



IF IKESET = 1 THEM CALL TGIMIT; 



CALL TGXHAC ; 



CALL GTP; 



CALL TGCOHP; 



CALL TGTRAW; 



CALL TGHZC; 


CALL TGSBC; 


CALL TGPHIC; 


A 


Figure 4.- TAEM guidance flow. 
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TGIMIT 

i 

i 

i 

i 

* 

I 

I PHASE = 1; 

I SR s RFTC / DTG; 
DMZCF = 0.0; 

QBLL = MXQBWT WEIGHT; 
DSBI = 0; 

PHILIM = PHILH1; 

DMZUL = DKZUC1; 


DMZLL = DHZIX1; 
OBARF = QBAR; 
QBD = 0; 

TGJEWD = 0 
I RESET = 0 

"I 

B 


Figure H 


- Continued. 
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B 

I 

TGXHAC 


1 

1 

IF Vf EIGHT > WT GS1 j 

THEM l 

i 

IGS = 

2; 

, 

i 

* 

: ; 

i 

i 

• 

i 



i 

* 

i 

i 

i 

i 

i I 

* 

■ • 

• 

t_ 

IGS s 

i 

i; I 

j 


IF GI_CHANGE = 1 THEN IGI = 2; 
ELSE Id = 1; 


XFTC = XA igi + HFTC igs / TGGS igs ; 
XALI = XA IG j + HALIjQg/ TGGSjjjg; 


XHAC = XFTC 


SHPLYK = -XHAC + DSHPLY; 
T 

i 

C 


Figure 4. - Continued. 
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TGCCHP 


ERPRED = RPR ED + XALI 
EOf = H + V V / (2G); 


\ 

IF UKPKED < EOW_SPT igs > 

/ 


THEN i 

! I EL = 2; 

I 


! ELSE 



EN s EN_C1 IGS J[EL ♦ (DRPRED - RN1 IGS ) EN_C2 IGS> IEL ; 


D 


Figure 4.- Continued 
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! THEM 


IF DR PR ED > PBRC igs !- 


HREF s ?BHC igs + PBCCjcs (DR PR ED - PBRC igs ) 


ELSE 


HREF = HALI igs - TGGS igs DRPRED 


\ THEN 

IF DRPRED > 0. > 

/ 


HREF = HREF 4 - DRPRED DRPRED (CUBIC_C3xcS * 
DRPRED CUBIC_CM igs ); 


\ THEN 

IF DRPRED > PBRCQjgs > 


QBREF = MIDVAL (QBRa IGS 4 - OBCIjqs (DRPRED - PBRCQjqs), 
QBRLL igs , CBRHL IGS ) 


ELSE 


OBREF = MIDVAL (QBRUL JGS 4. QBC2 ICS DRPRED. QBRLLr^. CBRULjnO 


IGS* wnui -IGS J 


HERROR = HREF - H 


IF DRPRED > PBRC igs 


DELRhG = HERROR/DHDRRF 


THEN 


DHDRRF = PBGC I0S 


ELSE 


DHDRRF = MIDVAL (-TGGS IGS + DRPRED (2 
CUBIC_C3 igs 4. 3 CUBIC_C4 IGS CRPRED),PBGC igs , - 

tggs igs ); 


CBARD = MIDVAL (COG (QBAR - CPARF) , - 
QBARDL, CBARDL); 

OBARF = QPARF + OBARD DIG; 

QBD = CDEOD QBD 4 - CQDG QBARD; 

QBERR = QBREF - QBARF ; 

EAS CMD = 17.1865 SORT (QBREF); 


Figure 4.- Continued. 



TGTRAM 


IF I PHASE = 3 


THEM 


IF (ABS(HERROR) < HERROR AND ABS(T) < 
TJERROR AMD ABS (GAMMA - GAHSG5 IGS ) < 
GAMHA_ERROR AMD ABS(QBERR) < QBERR0R1 
AMD H > HREF1 ) OR (ABS(HEFROR) < (H 
DELPHI - DELJI2) AND ABS(Y) < (H Y RANGE 1 - 
YJIAMGE2) AND ABS(GAMNA - GAKSGS ICS ) < ( 

H GAMMA JCOEF1 - GAMHACOEF2) AND ABS( 

QBERR) < QBERROR2 AND H < H REF1 ) OR 
H < H REF2 



/ 

/ 

/ 


TG DID x 1 


Figure 4.- Continued. 
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ELSE 


IF R PR ED < SHPLYK 


THEN 


I PHASE r 3; ! 

PHIO s PHIC; » 
PHILIH = PHILH3; ! 
DNZUL = DNZUC2; ! 
DNZLL = DIG LC2; I 


ELSE 


1 1 

I 1 

i : • 

J DO CASE IPHASE ♦ 1; ! 

i i 

i i 

: 26 ; 
* i i 

i « * 


S-TURN 




RETURN 


HA 


Figure *».- Continued. 
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25 


1 

1 

I 

I 

i*" 

* 

THEM 1 

I PHASE = 1 ; S 

IF EOT < EH *, 

1 

1 


t 

t 

1 

1 

1 

PHILB 1 s FHILM 1 ; | 


Figure - Continued. 
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27 


ES = ESIjqs + (DRP8ED - HHD<ST IGS ) EDRSjQg; 


I 

t 

I 


1 


THEM f 

• 

I 


IF ECW > ES AMD DR PR ED > BHIHST IGS { 

I 


• 

— "■ ~ i 

• 

I PHASE = 0; ! 

• 


1 

t 

1 


i 

t 

i 

* 

PHILD4 s PHILMO; J 

i 


1 

1 

1 


i 

• 

« 

i 

i 

S s ISQt ; » 

* 


1 

1 

1 

1 

I 


i 

• 

i 

• 

i 

SPSI = S PSD; 5 

• 

i 


1 

f 

f 

t 


i 

c 

• 

1 

■ a 

1 

THEM 

t 

1 

t 

1 


a 

i 

i 

IF SPSI < 0 OR SPSI > t35 ! 

1 

a 

t ' 

• 

i 

i 

» 

t 

1 

1 

1 



a 

i 

i 

• 

• 

» 

i 

i l 

i 

J S s 
• 

i EMEP = EMEP_C1 IGSi iel + (DRPRED - R»I1 IGS ) 

1 

t 

1 


i 

t 

I 

i 


1 

1 

•* 

» 

i THEM ! 

• 

• 

! EHEP_C2 IGSIEL ; 

1 

1 

I 

i 

i 

i 

i 

IF MEP s 1 J ! S = 1 

! J 

a 

• 

i 

1 

1 

1 

1 

1 

1 






Figure 4.- Continued* 
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&mas 


IF RCIR < P2TRHC1 RTURH 


THEM 


RT = (IDOT JCIR + TDOT TCBO / V; 


IF (RCIR < P2TRMC2 RTURM OR RT < RTBIAS) 


i 

i 


s 

s 


THEM 


S 

IPIUSL. 2 2; J 

philjw * rail* 2; ; 
t 


Figure *t.- Continued. 
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TGHZC 

I 

I 

I 

I 

GDH = MIDVAL(GDHC - GDHS H, GDHLL, GDHUL); 
HDREF r VH DHDRRF; 

HDERR = HDREF - HDOT; 

DMZC = DHZCG GDH (HDERR * HDREQG GDH 
HERROR); 

QBMHHZ = QBLL / AHA XI (COSPHI, CFHIM); 


QBHXHZ = QBHX1; 


THEN 


QBNZUL = -(QBG1 (QBWNNZ - QBARF) 

QBG2; 

QBHZLL z -(QBG1 (QBWXNZ - QBARF) - QBD) 
QBG2; 


IF HACH > QBH 1 


QBHXHZ z MIDVAL(QBHX2 + QBHXS (HACH - 
QBK2), QBHX2, QBHX3 ) ; 


- QBD) 


Figure 4. - Continued 



i ELSE 


EHAX = EK + EDELHZ IGS HIDVAL( (DS PRED - 


I 


R * 1 ICS > / DEL_B_EMiX IGS , EDELC1, EDELC2); 

EM IK = EM - EDELNZ 1GS ; 

EC&fMZUL s (GEUL GDH (EMAX - EOW) + HDERR) GDH 
GEHDUL; 

EGWMZLL = (GELL GDH (EMIM - ECM) + HDERR) GDH 
GEHDLL; 

DNZCD = (DMZC - DNZCF) CQG; 

DMZCF = MID VAL (DMZCF + DMZCD DTG, SCWKZLi., 
ECWMZUL); 

NZC = MIDVAL(DKZCF, CJ3NZLL, CSKZUL); 


IF I PHASE = 0 



GCOMT = (1. - QBARF / QBMXNZ) TAS / (( 
TAS 4- VCO) SECTH COS PHI); 

NZC = NZC - GCOMT; 


NZC = MIDVAL(NZC, DNZLL, DNZUL); 


Figure M.- Continued. 
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TG5BC: 



THEM 

} DSBC iff = DSBSUP; 


Figure 4. - Continued 
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scntas 


THEN ; : 

J 35 s 



8CFK25 


35 


DS8I = HIDVAL(D6BI + GS8I QBEftB DTG, - 
DSBIL, DSB3L ) ; 


Figure 4.- Continued. 
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TGPHIC 


I 

I 

I 

J 

if mach > mm 


THEM [ 

! PHILIHIT = FHILMSfJP; 


ELSE i 

1 PHILIHIT s FHIL3H; 

t 


DO CASE I PH AS: + 1; 


1 


t 

1 

1 

PH3E - S PHILIHIT; { 

J ' 

1 

I 

I 

i 2 

! i 

* i 

i 

i 

t 

t 

i 

i 

t 

i PHIC = GPHI DPSAC; J 

! i 

■ t 

! 3 


I 

1 

1 

RDOT = -C3CIR XDOT + TCIR TOOT) / BCIB; 

1 

1 

1 

I 

BERRC = HIDVAL(GR (RCIR - RTUHM) , -RERRLH, 

1 

1 

t 

1 

RERRLH); 

1 

1 

I 

1 

PHIC = TSGH (PHIP2C + RERRC + GRDOT RDOT); 

t 

1 

1 

1 


1 * 



« 

1 

TERRC s MIDVAL(-GT T, -TERRLH. TERRLH); 



PHIC = TERRC - GTDOT TOOT; 



1 ' 

1 

1 

1 

I 

1 


Figure Continued. 
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i 

IF ISI > 0 ? 


the* 


m& s (FHIE - mio) /ISB 
ISE 3 ISB - 15 
mi£ = PtilO * DFHX | 

PHIO = FfilC ; 


PHIC AT s HIDVALCPHIC, -PHILIMIT, FHILIHIT); 


» 

» 

i 

t 

§ 


Figure 4.- Continued. 
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mum 


GTP 


IF IPHASE = 3 J 

t 

r 

! 

I THEM 


RPR ED = SQRT(X I + TT); 


IF I PHASE < 2 


IF MEP = 1 


TSGM c SIGK(T); 


TSGM = -1, 
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THEN ; : 

! 43 : 

I I 

I I 

ELSE ! 

! RTAM = 0. ; 

t 

I 

I 

PSC r ARCTAM2(YCIR, XCIR); j 

t 

PST = (PSC - YSGH ARCTAM2(RTURH, RTAM)) ] 

I 

I 

RTD; ; 

i 

PST = RES 180 (PST); ! 

DPSAC = RES 180 (PST - PSD); \ 

« 

PHAVG = MIDVAL(PHAVGC - PHAtfcS MACH, | 

i 

PHAVGLL, FHAVGUL); J 

t 

RTAC = V VH / (G TANCPHAVG DTR ) ) ; ! 

i 

i 

ARCAC = RTAC ABS(DPSAC) DTR ; • 

i 

A = RTAC (1. - COS(DPSAC) ) ; | 

i 

B = RTAN - RTAC ABS (SIN (DPSAC) ) ; • 

• 

RC = SQRT(A A + B E) ; ! 

t 

PSHA = ABS(PST); J 


i 

J RPR ED = RCIR PSHA DTR - XHAC; 
* 

i 

ELSE | J 

! IF RCIR > RTURM ! 


\ 

IF XCIR < 0. AND ABS(Y) < 2. RTURN \ THEN 

OR ( YSGN SIGN (Y) < 0.) / 

/ 

v 

! Figure 4.- Continued. 
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I 

I 


ARCHA = RTURM PSHA DTR; 

RPRED s ARCAC * RC * ARCHA - BIAC; 


Figure 4.- Continued. 
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43 


RTAN = SORT (RCIR RCIR - RTURM RT' ); 


Figure 4.- Continued. 
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44 


PSHA = 360. - PSHA; 


Figure 4.- Concluded. 
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°.OFM? n , 


EH TER 

t 

i 

i PHIR * PHI*DTR i 

i SIHPHI * SIN(PHIR) | 

i COS PHI * COS (PH IP) ! 

I TANPHI « SINPHI/COSFHI | 

! COSTO * COS(THETA*DTR) \ 

\ ALFR « ALPHA#DTR ! 

! SINALF ■ SIH(ALFR) ! 

i COSALF • COS (ALFR) ! 

IlOil « MQD(LOOP, ?) | 

~ ] ■■■ - ■ — ■ 

| ROLL CHANNEL 

I 

1 

i GPBANK s MIDVAL (GPS*MACri + GPI,GPLL,GPUL ! 

! PCLIM = H IDVAL (PCS *KACH + PCI, POLL, PC UL) | 

! CALL SMOOTH (1,1 1 OH, LOOP, DTG, DT2, PHIC .BAHKSM) ! 

! BANKER = BANKSM-PHI ! 

! PC s MIDVAL( BANKER *GPrANK, -PCLIM, PCLIM) ! 

I BRATE(I) s PC* DTP ! 

■ ; ■■ 

i 

! PITCH CHANNEL 

i 

I XIN 7 R* MIDVAL (TANPHI, -T PLIM , TPLIM ) J 

1 CALL FILTER (1, LOOP, GF 11, GF IP, GF 13, DTP, XIN, RTANP) ! 

| CALL SMOOTH (2, II Oil, LOOP, DIG, OT2, NZC, NZCSM) ! 

! DNZCMP = ~CCSTH/MJDVAL(COSPHI,CPMIN, 1 . ) I 

! NZERR = NZ-( 1 , + VCO/TAS )*NZCSM + DNZCMP | 

— — • — 1 1 • 

i 

I 

I XIN = NZERR *GQN ! 

! CALL FILTER (2 , LOOP, GF21 , GF22, GF23, DTP, XIN, QC) ! 

! XIN = QC-RTANP ! 

! CALL FILTER (3, LOOP, GF31 , GF32, GF33, DT2, XIN , PCSL) j 

1 ALPMIN = MIDVAL ( AMNS*MACH + AMNI , AMNLL, AMNUL) ! 

""“I '■ r " ' ' ■"■■■'' ■ 1 “ ' " 

i 

i 

! IF \ ! ALPMAX = MIDVALCAMXS 1*(MACH - AMXM1) { 

1 MACH < AMXM2/ ! + AMXLL, AMXML, AMXLL) \ 

i ~~ — — ~ 

i 

i 

! ’ \ 

! IF \ ! ALPMAX = MIDVAL (AMXS2*( MACH - AMXM2) ! 

! MACH > AMXM2 / | + AMXML, AMXML, AMXUL ) I 

i “ / *“ 


A 


Figure 5.- TDAP detailed formulation. 
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A 

1 

i msrrgsp -? kmaimm 

I DQHIA * GQA*(ALPHA-ALPMIN) I 
I DQLOH « GQNL*<NZ-NZMAX) I 
I DQHXM « GQNL * ( HZ -NZH IK } i 

I DQLO * AMAX1 (DQLOA, DQLON) I 
I DQHI » AJ4XN1 (DQHIA, DQHIN) I 

fl?' \ THEM i BCSL * M IDVAL < BCSL , DQLO, DQlifV \ 

] DQLO < DQHI / " 

I 

i ^mmrrr^smr i 

~i 

i YAW CHANNEL 

t 

I 

| GRAY s M IDVAL (GYS*MACH + GYl , GYLL , GYUL ) i 

I DRPRM s COSTH*SINPHI*GRPHI/TAS j 

I RSTAB s DR PRM*COSALF-P*SINALF ! 

! DPC a GBDOT*RSTAB ! 

I CALL FILTER (4 , LOOP, GF ill, GFH2.GF M3, DY2, NY, NYF) | 

I DRRC * NYF*GRAY + DPC ! 

! BRATEC3) * ~DRRC*DTR ! 

‘ r 1 " " ' “ 11 " 1 ' " 1 “ 

i 

I ATTITUDE INTEGRATOR 

i 

i _______ 

! CALL EIGEN (BRATE, DT2, CBI, CBI) | 

I P a PC i 

i Q a -BCSL ! 

! R * -DRRC I 

""T 1 — — * 1 

i 

I SPEEDBRAKE CHANNEL 

i 

i 

I HR a M IDVAL (HMSB/HMS, 0. , 1 V) | 

! RMAX = SBRLC*SQRT ( 1 . -HR ) ! 

! RMAX a AM INI (SBRLO.RMAX) ! 

| RMIN s -SBRLC ] 

~i 

,i _________ 

! \ 

i IF \ THEN ! RMIN = -SBRS ! 

| DELB < SBSOFT / 

! “ / 
i 
i 

! DDOTRC = ( DSBC -DELB) /DT 2 ! 

IDDOTRC = M IDVAL (DDOTRC, RMIN, RMAX) ! 

1DELB = DELB + DDOTRC*DT2 ! 

i" 

i 

B 


Figure 5.- Continued. 
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B 

! 

» 

\ 



DELB * HIDVAL(DELB t SBHIH f DSBLIM) 
LOOP * LOOP + 1 



I 

I 

t 


END 


Figure 5#- Concluded. 
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ENTER 


I 

INPUT; NF, LOOP,C 1 ,C2, A, OT2, XIN 
OUTPUT; XOUT 


i 

i 


i \ 

I IF \ 

I LOOP • 0 / 

I / 

I 


I 

I 


THEN 


GXKNF) * Cl ! 

GX2(NF) + <1.-EXP(-A»DT2))*C2/A-C1 
CX3CNF) • -EXP(-A*DT2) ! 

XI (NF) * XIN ! 

XO(NF) * (GXKNF) + GX2(NF) )/( 1 . + GX3(NF) )*XIN { 

» 

. . I 


| XOUT « GX1 (NF) *XIN + GX2 (NF)*XI (NF)-GX3 (NF)*XO (NF) 
I XI (NF) s XIN 
! XO (NF) s XOUT 

i 

i 

~i — — — 

i 


RETURN 


Figure 6.- FILTER detailed formulation. 
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ENTER 

I 

1 

INPUT: 

OUTPUT 

I 

I 

I 


NS, IS, LOOP, THAJOR, TMINOR, XIN 
: XDUT 


IF 

LOOP 


\ 

\ 

r 0 / 

/ 


THEN ! XINP(NS) x XIN 


IF 

IS = 


\ 

\ 

0 / 
/ 


RETURN 


THEN 


ELSE 


XOUT « XI M 

XI NCR (NS ) s (XIN-XIHP(NS))»'I INOR/T'iA JOR 
XI NP (NS ) s XIN 



Figure 7.- SMOOTH detailed formulation. 
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EM TER 

i 

IMPUT: W,BT2,CB:i 

OUTPUT* CBI2 

i 


HAGW * VMAG(WE) 


■ y - ■ 

l 

-U 


IF 

MAGW 


THEN 


> 0,002 / 

/ 


r 


K 1 * ( 1 , ~C0S(MA'tf*DT2) )/MAGW**2 
K2 = SIN(MACW*rr.?)/MAGV/ 


ELSE 


K 1 s DT2*DT2/2. ! 
K2 = CTT 2 , 


DEC 1. 1 ) 
DEC 1,2) 
DE ( 1 , 3 ) 
DE (2, 1) 
DEC 2, 2) 
DEC2, 3) 
DE ( 3, 1 ) 
DEC3.2) 
DEC3,3) 


1,-(K1*(WE(3)**2) 
(K 1 *WE(2 )*WE( 1 ) )- 
(K1*WE(3)*WE(1 )) 
(K1 #WE(1 )*WE(2 ) ) 
1.-(K1*(WE(3)**2) 
(K1*WEC3)*WE(2))- 
( K 1 *WE ( 1 )*WE (3 ) )- 
(K 1 *WE (2 )*WE (3 ) ) 
1.-(K1*(WE(2)**2) 


)- (K1 *CWE (2 )**2 )) 
(K2*WE(3 ) ) 

+ (K2*WE(2)) 

+ (K2*WE(3)) 

)~ (K1 *(WE ( 1 ) **2 ) ) 
(K2*WE(1 )) 
(K2*WE(2)) 

+ (K2*WE(1 )) 
)~(K1*(WE(1 )**2)) 


I 

I 


CALL MAMPCCBI1, DE.CBI2) 

CALL CROSS CCBI2C 1,2) ,CBI2C 1,3) ,CBI2C 1, 1 ) ) 
CALL CROSS (CBI2C1 , 3),CBI2(1 , 1 ) ,CBI2( 1,2)) 


\ 

DO FOR I = 1 , 3 > 
/ 

( 

i 


CAl-l UNIT ( CBI2 ( 1,1), CBI2 (1,1) 


RETURN 


Figure 8.- EIGEN detailed formulation. 
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APPENDIX 

GUIDANCE AND FLIGHT CONTROL INPUT 
PARAMETER CALCULATIONS 


A-1 
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APPENDIX 

Guidance and Plight Control Input Parameter Caloulationa 

This appendix is included to clarify implementation requirements and oontains 
supplementary information which deeoribes the derivation of the TAEM guidance 
and flight control input parameters which are listed as GCQMP source inputs. In 
the flight software operating environment, these parameters are supplied by the 
attitude processor, air data, and navigation subsystems. Since these subsystems 
are not inoluded in the MCC simulation, their funotion must be approximated by 
calculation of the necessary parameters from the error-free vehiole state. Defi- 
nitions of these parameters, Which art derived in the following seotion, are 
contained in the input tablea of this document, 

a. State vector parameters - Sines ths Earth relative coordinates of the land- 
ing runway are constant, ths rotation matrix from Greenwich to runway co- 
ordinates, and the position vector of the runway threshold may be calcu- 
lated only on the initial pass and saved for use throughout the trajectory. 
This matrix is a Z, Y, Z Euler rotation sequence through the angles X RW , 

+ 90), and respectively. 


CER 


~-s<i>cXciP - sAcil; 

stbcXsMJ - sAc^ 
,-C<|>CA 


-S<t>SACi|J + CAC'P C\|JC4> 
S^SASty + CACty -CttSvl' 
~Cc>SA -S4) 


RREjj 


RREy 


RRE Z 


REQ 

■ ■ ^ + 

•%/ C4>^ + (1 - e ) 2 Sd> 2 

REQ 

.V^ 2 + (1 I - e) 2 Sd> 2 

REQ (1 - e) 2 

— r T | --» ■ - - - — ^ ^ 

+ (1 - e) 2 S(i} 2 


h RW 


h RW 


hRW 


Cd>CA 


CcpSX 


s«t> 


where 


s4> 

= sin 

SA = sin A r ^ 

= 

sin ^ Rv j 

Cep 

s COS $ RW 

CA = cos A rw 

Cty = 

cos Mj r ^/ 


The position and velocity vectors in Earth relative Greenwich coordinates are 
calculated using the time dependent inertial to Greenwich transformation matrix. 


.vi-.unjiN’t; iy.ci.v 
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If the inertial coordinate system used is based on the truu date of verno l 
equinox and Earth polar axis, this matrix is simply a Z-ax's rotation through 
the Greenwieh hour angle, and the Earth angular rate vector is coincident with 
the Z-axie. 


RE * 


» 


VB * 


CIE HI 

: s __ 

CIE {VI - w e x HI) 


The position and velocity vectors in runway coordinates ar‘'j 


RRW = 


CER 


(RE - RRE) 


VRW = 


CER 


VE 


To construct the topodetic coordinate system, the longitude and geodetic 
latitude of the vehicle must be obtained. This is done by first determi:: : 
the quantity A by iterating three times the equation 


A = 


+ REQ 


1 - (1 - e) 2 

\/RXY 2 /A 2 + (1 - e) 2 RE z 2 


where 


RXY = /RE X 2 + RE y 2 


The starting value used for A is 1/(1 - e) 2 1.0067, which is its true value 

for h = 0. Three iterations will guarantee about 7 to 8 digit accuracy. 

The geodetic latitude and altitude are given by 


<}> D = ATAN 


RE Z 

. (RXY/A). 


»D 


1 - A(1 - e) 2 , s — 

VRXY 2 /A 2 + RE Z 2 

1 - (1 - e) 2 
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Longitude it easily obtained from 
X s ATAM2(RE y , RE X ) 

The Greenwich to topodetio tranafortaation matrix is a Z, Y, Euler rotation 
sequence through the angles X, and + 90), respectively. 

^0 ooa X -ain sin X cos 4 >d * 

X oos X 0 

oos X -cos ♦p sin X -sin $d 

related parameters can now be obtained. 


VH *-\/VT x 2 + VT y 2 

v a |vt| 

GAMMA = ASIN (HDOT/V) ' RTD 

H * h D - h RW 
X + RRW X 

X = RRWy 
XDOT = VRW X 
YDOT a VRW y 

PSD = ATAN2( VRWy , VRW X )'RTD 


CET 


•sin 
•aln 
L-eos 


The topodetio velocity 

H 


VT a 
HDOT 


VE 


-VT, 


b. Attitude parameters - The pitch and roll angles are obtained from the body 
to topodetio transformation matrix. 


CBT . 
THETA 


[cet] [cie] [cbi] 

-ASIN(CBT( 3j1 )) - RTD 
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PHI « ATAN2(CBT(3 } 2) j CBT ( 3 , 3 > 3 'RTD 
COSPHI * COS(PHl/PTD> 

SECTH « 1 ,/CQS(TKETA/RTD) 


The Shuttle flight control sensors include body axis accelerometers and rite 
gyros which provide translational acceleration and rotational rate measurements 
to the digital autopilot. Por this simulation, these will be assumed to be per- 
fect systems providing true quantities. Since body rates are calculated in 
TDAP, they will be available as measured values for the following pass. The 
body axis aoosleroaster measurements are obtained from the true aerodynamic 
f cross . 


NY s FABy /WEIGHT 
NZ s -FAB Z /WEIGHT 


c. Air data parameters - The Shuttle avionics includes an air data subsystem 
which derives air relative free stream parameters (mach number, true 
airspeed, dynamic pressure, and anglo of attack) from the pilot-static probe 
pressure measurements. Because of thermal constraints, the air data probes 
oannot be deployed until approximately mach 3.5. Air data parameters will 
be supplied to the TAEM guidance and flight control below a predetermined 
mach number (MACHAD) which will be mission dependent. This mach number will 
be 2.5 for the STS-1 mission. The parameters supplied by the air data sys- 
tem are required continuously by the TAEM guidance and fight control and, 
therefore, must be obtained elsewhere when acceptance of air data is undesir- 
able. In the flight software, this is done Dy deriving the required 
quantities from the navigated state vector and attitude data, assuming calm 
wind conditions and standard atmosphere. The inaccuracies which result in 
the presence of winds and nonstandard atmosphere have an effect on guidance 
performance and will be included in the MCC TAEM simulation by assuming per- 
fect navigation and air data systems, with the changeover occurring at the 
proper time. This changeover should be irreversible so that air data will 
be used continuously after its first acceptance. 

For the region prior to acceptance of air data (MACH > MACHAD), the velocity 
of sound is assumed constant (1000 fps) and standard air density is 
approximated by 


(-h D /20600. ) 

p 3 = . 00413579 e (h D > 35000.) 

C-h D / 30550.) 

P 3 s .002376B9 e (h D < 35000.) 
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To oaloulate angle of attack in this region, the Earth relative velocity 
vector aust be obtained in body coordinates. 

HI • [cbt] t W 


The required quantities for the no air data region can now be defined. 
TAS i |v5f| 

ALPHA • ATAN (VBE Z /VBE X ) 'RTD 
MACH * TAS/ 1000. 

QBAR * 0 ,5P S TAS 2 


After the acceptance of air data (MACH < MACHAD) , the corresponding true 

air relative quantities are to be sent to the TAEM guidance and flight control. 

vbw s [cbt]t (VT - VW) 

TAS s | VBW ( 

ALPHA = ATAN (VBW Z /VBW X ) • RTD 
MACH = TAS/ VS 
QBAR = 0.5 PA TAS 2 
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CBI 

CBT 

CER 

CEX 

CIS 

0 

FAB 

hp 

h RW 

MACHAD 

RE 

REQ 

RI 

RRE 

RRW 

RTD 

VBE 

VBVf 

VE 

vT 

VRW 

VS 

VT 

vw 

A 


SYMBOLS 

body to inertial coordinate transformation matrix 

body to topodetic coordinate transformation matrix 

Greenwioh to runway coordinate transformation matrix 

Greenwich to topodetio coordinate transformation matrix 

inertial to Greenwich coordinate transformation matrix 

ellipticity (flattening) of Earth ellipsoid model 

total aerodynamic foroe vector (body coordinate) 

geodetic altitude of vehicle center of mass 

geodetic altitude of runway threshold 

true mach number to initiate acceptance of air data 

position vector of vehicle center of mass (Greenwich coordinate) 

equatorial radius of Earth ellopsoid model 

position vector of vehicle center of mass (inertial coordinate) 
position vector of runway threshold (Greenwich coordinate) 
position vector of vehicle center of mass (runway coordinate) 
radians to degrees conversion constant 
vehicle Earth relative velocity vector (body coordinate) 
vehicle air relative velocity vector (body coordinate) 
velocity vector of vehicle center of mass (Greenwich coordinate) 

velocity vector of vehicle center of mass '’inertial coordinate) 

velocity vector of vehicle center of mass (runway coordinate) 

true atmospheric velocity of sound 

velocity vector of vehicle center of mass (topodetic coordinate) 
Earth relative velocity vector of air mass (topodetic coordinate) 
longitude of vehicle center of mass 
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longitude of runway threshold 
Pa true atmospherio air density 

P 8 standard air density estimate 

♦d geodetio latitude of vehicle center of mass 

$RW geodetic latitude of runway threshold 

runway azimuth with respect to true north 
Earth angular rate vector (inertial coordinate) 
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